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25 ABSTRACT 


26 Water availability plays an important role in the socio-economic development of a region. It is 

27 however, subject to the influence of large-scale circulation indices, resulting in periodic excesses 

28 and deficits. An assessment of the degree of correlation between climate indices and water 

29 availability, and the quantification of changes with respect to major climate events is important for 

30 long-tenn water resources planning and management, especially in transboundary basins as it can 

31 help in conflict avoidance. In this study we first establish the correlation of the Pacific Decadal 

32 Oscillation (PDO) and El Nino-Southern Oscillation (ENSO) with gauged precipitation in the Rio 

33 Grande basin, and quantify the changes in water availability using runoff generated from the Noah 

34 land surface model. Both spatial and temporal variations are noted, with winter and spring being 

35 most influenced by conditions in the Pacific Ocean. Negative correlation is observed at the 

36 headwaters and positive correlation across the rest of the basin. The influence of individual ENSO 

37 events, classified using four different criteria, is also examined. El Ninos (La Ninas) generally cause 

38 an increase (decrease) in runoff, but the pattern is not consistent; percentage change in water 

39 availability varies across events. Further, positive PDO enhances the effect of El Nino and dampens 

40 that of La Nina, but during neutral/transitioning PDO, La Nina dominates meteorological conditions. 

41 Long El Ninos have more influence on water availability than short duration high intensity events. 

42 We also note that the percentage increase during El Ninos significantly offsets the drought-causing 

43 effect of La Ninas. 

44 Keywords: Climate variability, PDO, ENSO, MEI, EMI, Noah LSM, Rio Grande, Water 

45 Availability 
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INTRODUCTION 


48 Large-scale climate patterns have a significant influence on local atmospheric and hydrologic 

49 variables, and consequently on water availability. Several studies have investigated the influence of 

50 climate variability using either a single index or a combination of indices on precipitation [e.g., 

51 McCabe and Dettinger, 1999; Piechota and Dracup, 1996; Ropelewski and Halpert, 1986; 

52 Woolhiser et al., 1993], streamflow [e.g., Barlow et al, 2001; Kahya and Dracup, 1993; Redmond 

53 and Koch, 1991], and drought [e.g., Ozger et al., 2009; Schoennagel et al., 2005] in the United 

54 States (US). In the southern US, the Pacific Decadal Oscillation (PDO) and the El Nino-Southern 

55 Oscillation (ENSO) have been found to be the two most dominant climate teleconnections 

56 influencing regional hydrological conditions. 

57 ENSO is a coupled ocean-atmosphere phenomenon related to sea surface temperature (SST) 

58 anomalies (SSTA) in the central and eastern equatorial Pacific and associated sea-level pressure 

59 difference known as the Southern Oscillation [Rasmusson and Carpenter, 1982]. It has a recurrence 

60 pattern of 3 to 6 years and every event normally lasts for about a year. El Nino events, the positive 

61 or wann phase of ENSO, are often, but not always, followed by La Nina events, also referred to as 

62 the negative or cold phase of ENSO. More recently, a new type of El Nino, occurring more 

63 frequently, with inter-annual variability, has been observed [Ashok et al., 2007]. It has been named 

64 El Nino Modoki (Japanese for “similar but different”). It occurs with a shift in the warming center 

65 from the eastern equatorial Pacific, which is the case with regular El Nino, to central equatorial 

66 Pacific, and both the eastern and western regions are flanked by anomalously cool temperatures, 

67 thus resulting in an SST gradient that generates a two-cell Walker Circulation in the troposphere 

68 with a wet region over the central Pacific. When coupled with other ongoing atmospheric 

69 disturbances, a dry rim arises around the wet central tropical Pacific. Given the similarities between 

canonical ENSO and this new occurrence, it is easy to confuse between their impacts. Partial 
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71 correlation and regression analyses suggest that they are distinct phenomena in both space and time 

72 and do not appear as an evolving phase of one or the other [ Weng et al., 2007]. 

73 PDO is a long-lived El Nino-like pattern of Pacific climate variability with a cycle of about 20 to 30 

74 years [Mantua and Hare, 2002]. PDO influences the hydrological cycle in the same way as ENSO, 

75 but with more pronounced influence in the extra tropics and secondary influence in the tropics. The 

76 similarities in the signature between ENSO and PDO have led to the hypothesis that the two 

77 teleconnections may be related, or PDO may be forced by ENSO [ Zhang et al, 1997]. Statistical 

78 analysis by Newman et al. [2003] showed that PDO is dependent on ENSO on all timescales. When 

79 the PDO is in its positive or wann phase, above nonnal SST is observed along the west coast of 

80 North America and below normal SST along the central and western North Pacific around 45°N. 

81 The Aleutian low strengthens and winter precipitation increases in the southern US [ Mantua et al., 

82 1997], 

83 Recently Kurtzman and Scanlon [2007] examined the impacts of ENSO and PDO on winter 

84 precipitation in 165 climate divisions in southern and central US and found a significant increase 

85 (decrease) with respect to El Nino (La Nina). The correlation with PDO was weaker, but when both 

86 indices were combined, it was noted that La Ninas occurring during the cold phases of PDO 

87 exhibited strong influence in central US and El Ninos occurring during the wann phases of PDO 

88 dominated southwest and southeast US. Redmond and Koch [1991] noted that if events in the 

89 Pacific Ocean are causally related to remote meteorological variables, it would be separated by a 

90 time lag. They found statistically significant correlations with lags ranging between 0 to 6 months 

91 between the Southern Oscillation Index (SOI) and precipitation in the western US. Kumar and 

92 Hoerling [2003] found that the maximum correlation between observed zonal mean tropical 200-mb 

93 heights and SST in the Pacific occurs with a lag of 1 to 3 months and this results in a lag of one 
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94 season between rainfall in the tropical Pacific and Nino 3.4 SSTA. The robustness of these results 

95 was confirmed using an atmospheric general circulation model (GCM). 

96 While these information, on the degree of association between hydrologic variables and climate 

97 patterns, are valuable, they are only of qualitative nature and their use in water management is 

98 limited. Water planning and management is driven by demand, priority, and availability. Demand is 

99 influenced by demographic and economic changes, while priority is an institutional variable defined 

100 by legal, social, and economic constraints. Availability, on the other hand, is a natural variable 

101 subjected to the whims of climate. 

102 In this paper we investigate the influence of large-scale climate indices, namely ENSO and PDO, on 

103 the water availability within the Rio Grande/Rio Bravo del Norte basin (RG). RG is a transboundary 

104 basin shared between three states in the US and straddles US and Mexico, two countries very 

105 dissimilar economically. It is a vital source of water for the region, but is already in a state of 

106 absolute water scarcity, with less than 500 m /person/day; the only transboundary basin in this 

107 category [Wolf 2002]. This region is also extremely vulnerable to droughts; records show that it 

108 suffers from both short-tenn and long-term droughts [Quiring and Goodrich, 2008]. Subjected to a 

109 burgeoning population, which will further increase the stress on water allocation, and climate 

110 change, which will likely result in a decrease in precipitation [IPCC, 2007], the potential for 

111 conflicts cannot be overlooked. It is therefore imperative to understand the mechanisms driving 

112 water availability and quantify any change for long-term sustainable water planning and 

113 management. 

114 The study is divided into two main sections. We first establish the correlation between large-scale 

115 circulation indices and gauged precipitation to explore the spatial and temporal influence of ENSO 

116 and PDO separately on the basin. The correlation structure of precipitation with two canonical 
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117 ENSO indices (Nino 3.4 and the Multi-ENSO index (MEI)) were compared to that of the El Nino 

118 Modoki index (EMI) to determine which index shows maximum correlation and is best suited for 

119 water management within the basin. 

120 Runoff is not linearly related to precipitation, but affected by natural processes and subjected to 

121 other meteorological variations such as temperature, evapotranspiration, wind speed, etc. which are 

122 also influenced by remote climate teleconnections. Streamflow is often used as a measure of surface 

123 water availability, but in large basins, stream gage records are not a realistic representation of actual 

124 flow as they are affected by dams, diversions, return flows, reduction in base flows by excessive 

125 groundwater pumping, and urbanization [Legates et al., 2005], thus obscuring climate influences. 

126 We therefore use a land surface model (LSM), but keep land-use-land-cover constant, to generate 

127 runoff as it incorporates all necessary factors in the process. The basin is divided into six sub- 

128 regions and the temporal variations in water availability with respect to climate indices are 

129 examined. 

130 All ENSO events are not similar and coincide with different phases of PDO (positive, negative, or 

131 transitioning from one phase to another). We compare and rank individual El Nino and La Nina 

132 events based on their durations, maximum (or minimum) SSTA recorded, and intensities - a new 

133 metric that we propose in this study. The percentage change in water availability in each sub-region, 

134 with respect to individual El Nino and La Nina events and coincident PDO phases is then examined. 

135 The paper is organized as follows. Section 2 discusses the methodology and the choice of LSM 

136 adopted for this study. A description of the study basin is given in section 3. Section 4 describes the 

137 datasets for precipitation and climate indices, and forcing and parameters for the LSM. The first part 

138 of section 5 discusses the hydroclimatology of the basin, the correlation between climate indices and 

139 gauged precipitation, and the differences between ENSO events. In the second part we discuss the 
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140 model output and validation, and investigate the lags and changes in water availability with respect 

141 to climate variability. The main conclusions drawn from the study are given in section 6. 

142 2 METHODOLOGY 

143 Pearson correlation is used to detennine the relationship between climate indices and gauged 

144 seasonal precipitation data. Kriging interpolation is employed to map the spatial variation of the 

145 correlation coefficient across the study area. Runoff, which is a proxy for surface water availability, 

146 is obtained from an LSM. The factors considered in the selection of the most appropriate LSM for 

147 the region is discussed. Continuous wavelet transfonn is used to investigate the temporal structure 

148 and influence of climate variability on water availability. 

149 2.1 Correlation 

150 The Pearson correlation coefficient, p xy , is a measure of linear association between two time series: 

151 x and y . The lag-correlation coefficient, p xy (k), is the cross-correlation for lag k between the time 

152 series. The range for p xy (k) is [-1, 1], with larger p xy implying greater ability of x to predict y 

153 [von Storch and Zwiers, 2003]. The correlation coefficient can be used as a statistical test of 

154 independence to help make inferences about the degree of association between variables. The null 

155 hypothesis is that the two time series are independent and identically distributed (iid) nonnal 

156 random variables [p xy = o). 

157 2.2 Land Surface Modeling 

158 LSMs compute terrestrial water, energy, momentum, and bio-geochemical exchange processes by 

159 solving the governing equations of the soil-vegetation-snowpack medium [Peters-Lidard et a/., 

160 2007]. A number of LSMs have been developed over the last 30 years and are constantly being 

refined as our understanding of the physics underlying earth system processes improves, and 
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computing capabilities increases. Four LSMs, namely Mosaic, Noah, the Variable Infiltration 
Capacity (VIC), and Sacramento models, were evaluated over the Continental US (CONUS) as part 
of the North American Land Data Assimilation System (NLDAS) project [Mitchell et ai, 2004]. 
Lohmann et al. [2004] evaluated these models for their performance in partitioning water balance 
terms (evapotranspiration, runoff, and storage change) across four different quadrants over CONUS, 
and their ability to reproduce streamflow at different timescales (daily, monthly, and annual) and 
noted that at the continental scale the results varied significantly in the wet eastern US but were 
generally in agreement over the drier western region. The Mosaic and Sacramento models 
underestimate runoffs, and VIC produced more runoff, while Noah’s predictions fell in between. In 
small- to medium-sized catchments, the models showed similar bias gradients in the east, increasing 
from north to south. VIC, for example, produced the right annual runoff in the northeast US but 
more runoff towards the south. Noah predicts less runoff in the northeast US, more in the south, and 
the right runoff in the middle. 

Noah was also found to exhibit the lowest regional bias. Lohmann et al. [2004] further speculated 
that the high runoff produced by the models in the southwestern US may be attributed to farming 
and irrigation which is not included within the NLDAS setup. Seasonal analysis showed that Noah 
produced correct runoff in a number of basins for the cold season with the same north-south annual 
bias. Both VIC and Noah produced soil moisture anomalies close to observed values. However, in 
the Little River Experimental Watershed, Sahoo et al. [2008] found that Noah produced higher soil 
moisture, which, as a result of the model physics governing partitioning, produced less surface and 
subsurface runoff. Nevertheless Lohmann et al. [2004], found that the Sacramento and Noah 
reproduced daily streamflow better, with Noah having the highest overall score based on the Nash- 
Sutcliff efficiency. The study also evaluated the model performance over nine large basins in the 
US. RG was not included, but the sizes of the basins examined are comparable. There was 


6 



186 disagreement between modeled and measured runoff in high and less regulated basins. In high 

187 regulated basis, which comprise a number of dams and reservoirs, smaller seasonal signals were 

188 observed, whereas in less regulated ones, seasonality was closely captured. This confirms that the 

189 models were actually effectively reproducing seasonal variations, which is being dampened by 

190 engineering infrastructures in the highly regulated basins, thus lending credence in using an LSM in 

191 modeling RG which has a large number of dams, diversions, and reservoirs. 

192 2.2.1 Noah LSM 

193 The community Noah LSM’s [ Chen et al., 1996; Koren et al., 1999] legacy extends into modeling 

194 efforts carried in the 1980s [Mahrt and Ek, 1984; Mahrt and Pan, 1984; Pan and Mahrt, 1987] and 

195 has been further refined in the 1990s under the GEWEX/GCIP/GAPP Program Office of 

196 NOAA/OGP, led by the National Centers for Environmental Prediction (NCEP) and benefitting 

197 from the collaboration of investigators from both public and private institutions. Noah has been a 

198 candidate in major off-line land surface experiments, such as the Project for Intercomparison of 

199 Land-surface Parameterization Schemes [PIPLS; Henderson-Sellers et al, 1996] and the Global Soil 

200 Wetness Project [GSWP; Dirmeyer et al., 1999] among others. It has been validated in both coupled 

201 and uncoupled modes [Mitchell, 2005]. 

202 Noah is a stand-alone 1-D column model that can simulate soil moisture (both liquid and frozen), 

203 soil temperature, skin temperature, snowpack depth, snowpack water equivalent, canopy water 

204 content, and water and energy flux tenns of the surface water and energy balance [Mitchell, 2005]. 

205 The model has a snow layer and a canopy layer. The soil profile extends to a depth of 2 m divided 

206 into four layers from the ground surface to the bottom: 0-0.1 m, 0. 1-0.4 m, 0.4-1 m, and 1-2 m. The 

207 root zone is limited to the upper 1 m of soil, and the lower 1 m layer acts as a reservoir with gravity 

208 drainage at the bottom. The snow layer simulates snow accumulation, sublimation, melting, and heat 

209 exchange at snow-atmosphere and snow-soil interfaces. Precipitation is deemed snow if the 
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210 temperature of the lowest atmospheric layer is below 0°C. The total evaporation, in the absence of 

211 snow, is the sum of direct evaporation from the topmost soil layer [Mahfouf and Noilhan, 1991], 

212 evaporation of precipitation intercepted by plant canopy, and transpiration from canopy of 

213 vegetation [Jacquemin and Noilhan, 1990; Noilhan and Planton, 1989]. Surface runoff is the excess 

214 after infiltration [ Schaake et al., 1996], A complete description of the model physics and order in 

215 which computations are carried out is available in Chen and Dudhia [2001]. 

216 Lohmann et al. [2004] point out that one notable difference between Noah and other LSMs 

217 considered in the NLDAS project is the earlier onset of runoff in snowmelt season when compared 

218 to other models and observed values. This may make Noah a less likely candidate for streamflow 

219 studies in snow dominated watersheds. In Noah, snow can either sublimate or melt as there is no 

220 horizontal transport. Sheffield et al. [2003] and Pan et al. [2003] evaluated the four LSMs 

221 considered in the NLDAS project for snow cover extent and simulated snow water equivalent. 

222 Systematic low biases were observed in the snow cover extent and snow water equivalent in the 

223 simulations for all four models and larger discrepancies were observed at higher elevations. Noah 

224 consistently underestimates the snow cover extent at all elevations. This under-prediction is partly 

225 explained by its higher snow water equivalent threshold for large snow cover values as compared to 

226 other models. Noah also tends to melt snow earlier, which is due to the low albedo values in each 

227 snow covered grid, which leads to higher available energy at the surface creating a positive feedback 

228 mechanism which enhances snowmelt and sublimation. 

229 Hogue et al. [2005] evaluated the transferability of calibrated parameters in Noah between two 

230 semi-arid sites in the southern US for evaluating model perfonnance under the different climatic 

231 conditions these regions are subjected to. They found that generally Noah accurately simulated 

232 sensible heat, ground heat, and ground temperature. However discrepancies were noted during brief 

periods of moist air influx responsible for monsoon or during El Nino. When compared against 
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234 other LSMs Noah reproduced streamflow with high accuracy, with the smallest bias in both 

235 evaporation and runoff with respect to observed annual water budget [Mitchell et al., 2004]. 

236 Thus, based on extensive evaluation and comparison and despite some of its limitations, Noah 

237 seems the best suited LSM for the hydrological modeling of RG for purposes of this study. This 

238 study also provides an opportunity to test the validity of Noah outside the southern US border, while 

239 still within NLDAS-2 domain, and may thus supplement previous findings. 

240 2.2.2 Land Information System 

241 Noah LSM was run within NASA Goddard Space Flight Center’s (GSFC) Land Information System 

242 (LIS; http://lis.gsfc.nasa.gov/). LIS, designed as a problem solving environment for hydrologic 

243 modeling applications, is an integrated high-performance land surface modeling and data 

244 assimilation framework [Kumar et al., 2006; Peters-Lidard et al., 2007] which evolved from the 

245 Global Land Data Assimilation System [GLDAS; Rodell et al, 2004] and the North America Land 

246 Data Assimilation System [NLDAS; Mitchell et al., 2004]. 

247 2.3 Continuous Wavelet Transform 

248 Wavelet transform decomposes a signal in terms of some elementary functions derived from a 

249 “mother wavelet” using a sliding window function whose radius increases in space (i.e., decreasing 

250 in frequency), allowing the low-frequency content of the signal to be resolved [Rivera et al., 2004], 

251 A number of mother wavelets are commonly used, and can be grouped into continuous and 

252 orthogonal wavelets; each having intrinsic advantages for specific applications. We chose the Morlet 

253 wavelet, which consists of a plane wave modified by a Gaussian envelope, as it is the most widely 

254 used continuous wavelet in geophysical applications. Its complex structure allows the detection of 

255 both time-dependent amplitude and phase for different frequencies in the time series [Lau and 

256 Weng, 1995]. 
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257 Kumar and Foufoula-Georgiou [1997] discussed the applications of wavelet to geophysical 

258 processes. Torrence and Compo [1998] reported on the application of wavelet, including a 

259 comparison to the widowed Fourier transform, to climate analysis using ENSO series. Labat et al. 

260 [2001] showed how wavelet can be appropriately employed in rainfall-runoff analysis. 

261 3 STUDY AREA 

262 RG is the fifth longest river in North America. It originates in the snow dominated San Juan range in 

263 the Rocky Mountains in southern Colorado, at an altitude of around 3,700 m amsl, and flows 

264 through arid/semi-arid plains in a south-eastward direction over a distance of approximately 3,100 

265 km before discharging into the Gulf of Mexico. The basin encompasses an area of 557,722 km" 

266 straddling three states in the US and five states in Mexico (Figure 1). The river catchment is narrow, 

267 with its length being considerably longer than its width, and has a dendritic drainage pattern. The 

268 watershed contains a number of endorheic sub-basins, such that only 468,374 km (242,994 km on 

269 the US side and 225,380 km“ on the Mexican side) actually contribute to flow in the river [Patino- 

270 Gomez, 2005]. 

271 Figure 1: The Rio Grande/Rio Bravo del Norte basin 

272 Winter precipitation and spring runoff sustain flow in the basin. The flow is impounded in a number 

273 of dams and regulated by major diversions. Elephant Butte Dam in southern New Mexico supports 

274 agricultural production in the region. The release from the dam is apportioned between the US and 

275 Mexico under the 1906 Convention for the Equitable Distribution of Waters of the Rio Grande 

276 [ International Boundary & Water Commission, 1906], From El Paso to Ojinaga/Presidio the river 

277 flows through one of its driest stretches [US Army Corps of Engineers, 2008]. At El Paso (station 45 

278 in Figure 2), for example, the mean annual rainfall is 219 mm, while the yearly pan evaporation is 
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279 around 1,500 mm. At Ojinaga/Presidio, RG is regenerated from flow from the Rio Conchos, which 

280 is one of its most important tributaries, originating in the Sierra Madre Occidental in northwestern 

281 Mexico at around 3,500 m amsl. Two international reservoirs, Amistad and Falcon, store and 

282 apportion the water between the US and Mexico. Without any dam or diversion along its course, the 

-5 

283 virgin flow of RG is estimated at above 100 m /s [ Revenga et al., 2003]. However, the over 

284 anthropogenization of the basin has constantly impacted flow, such that the actual mean historical 

285 flow is 37 m /s, but in recent years the flow in the river has reduced considerably, and on several 

286 occasions the river failed to reach the sea. 

287 Land cover in RG is mainly desert shrubland and grassland, covering about 81%, while irrigated 

288 agriculture constitutes only 2.6% of the basin, and urban and industrial area covers 6% of the basin 

289 [Revenga et al., 2003]. 

290 Given its size, the varied climatology it is subjected to, and the major dams and diversions 

291 partitioning the system, RG cannot be studied as one watershed, nor can it be divided into sub- 

292 basins, as some tributaries, like the Pecos, runs along the main stem, crossing several latitudes, thus 

293 subjected to different climate teleconnection influences. Higher snowfall at the headwaters for 

294 example does not necessarily result in higher flow into the international reservoirs or at the mouth. 

295 Moreover, climate divisions from one US state do not align with that from another state. The basin 

296 is therefore divided into six sub-regions (Figure 2) by considering the latter constraints and other 

297 geomorphological features in the system. 
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298 4 DATA 

299 4.1 Precipitation 

300 4.1.1 United States 

301 The National Climatic Data Center (NCDC) has an extensive archive of publicly available weather 

302 data from NOAA’s Cooperative Observer Program (COOP) stations in the US. The COOP data 

303 contains gaps and only a few stations have continuous records. NCDC also houses the United States 

304 Historical Climatology Network (USHCN) version 2 dataset which is a designated subset of the 

305 COOP network [ Karl et al., 1990]. The dataset undergoes extensive quality control including 

306 adjustment for any time-of-observation bias. Only 27 stations out of the 1,221 stations in the 

307 USHCN are serially complete [ McRoberts and Nielsen- Gammon, 2011], while missing data in the 

308 others have been filled using a weighted average of values from highly correlated neighboring 

309 stations. However, the density of the USHCN network is not adequate within RG for purposes of 

310 this study (see Figure 1 of USHCN Version 2 Serial Monthly Dataset, available from 

311 http://www.ncdc.noaa.gov/oa/climate/research/ushcn/ which shows the distribution of both COOP 

312 and HCN sites). 

313 McRoberts and Nielsen-Gammon [2011] proposed the full network estimated precipitation (FNEP, 

314 available from http://atmo.tamu.edu/osc/fnep) which utilizes as many COOP observations from the 

315 network of more than 24,000 stations and an inverse distance weighting scheme, using stations that 

316 have at least 10 years of overlap data and highest correlation, to fill missing data and extend the 

317 record, thus creating a continuous series from 1895 to present. There are a total of 332 FNEP 

318 stations within the US portion of RG, but only those that have a sufficiently long record of 

319 observational data, while ensuring adequate spatial coverage, were chosen, hence taking advantage 

320 of the filled gaps. The time period considered in the analysis is January 1935 to December 2008, 
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321 thus giving 75 years of data. Figure 2 shows the spatial distribution of the 63 stations (4 in Colorado, 

322 40 in New Mexico, and 19 in Texas) selected. 

323 4.1.2 Mexico 

324 Historical monthly precipitation data for the Mexican section of the basin was obtained from the 

325 Servicio Meteorologico Nacional (SMN), Comision Nacional del Agua (CONAGUA), Mexico. 

326 SMN is the state entity responsible for the observation, recording, interpretation and dissemination 

327 of weather infonnation in Mexico. The most updated observational dataset obtained contained data 

328 up to December 2006 only, and any record beyond that date is not yet publicly available. For the Rio 

329 Conchos sub-basin the latest available data is only up to December 2003. 

330 Of all the stations within the Mexican section of the basin, excluding the Rio Conchos, 52 are 

331 operational and have records that extend up to 2006. However, most of these stations are very recent 

332 and only a few have records of at least 50 years, but have several years of missing data. A careful 

333 selection of viable stations, while ensuring maximum possible coverage, limited the number of 

334 stations that can be used in the study to only 12. The time period extended from January 1954 to 

335 December 2006. Records show that 9 stations are operational in the Rio Conchos, but only 4 had 

336 sufficiently long records, extending from 1964 to 2003, for this study. Data gaps were filled from 

337 neighboring stations that had sufficient amounts of overlapping data with the target stations. Any 

338 data missing from both the target and neighboring stations were filled with the long term monthly 

339 mean. The spatial distribution of the selected stations is shown in Figure 2. 

340 Figure 2: Location of precipitation stations used in this study 
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341 4.2 Climate Indices 

342 4.2.1 Pacific Decadal Oscillation (PDO) 

343 Monthly PDO indices for the period 1935 to 2008 were obtained from the Joint Institute for the 

344 Study of the Atmosphere and Ocean (JISAO, http://jisao.washington.edu/pdo/PDO.latest). The PDO 

345 index is defined as the leading principal component of the North Pacific, poleward of 20°N, mean 

346 monthly SSTA [Mantua et al., 1997]. The data is not influenced by global wanning trends as 

347 monthly mean global average SSTA have been removed. A plot of the monthly indices, with a 

348 centered 13-month moving average to highlight multi-decadal frequency is given in Figure 3. 

349 Positive (negative) values indicate wann (cold) phases of PDO. Between 1935 and 2008, 51% of 

350 warm months and 49% of cold months were recorded. 

351 Independent studies have shown two full PDO cycles in the last century [Mantua and Hare, 2002; 

352 Minobe, 1997], cool regimes that lasted from 1890 to 1924 and from 1947 to 1976, and warm 

353 phases that lasted from 1925 to 1946 and from 1977 to 1998. From 1998, PDO has been in a cold 

354 phase until 2002 and in a wann phase from 2002 to 2007. Occasional shifts, within the 20-30 year 

355 cycle, from cool (warm) to warm (cool) are visible in the record. 

356 4.2.2 Nino 3.4 

357 There are four Nino regions along the equatorial Pacific, chosen in the early 1980’s to describe and 

358 monitor SST. The wanning across this region is not unifonn and no single region can capture the 

359 whole ENSO phenomenon. Barnston et al. [1997] proposed the Nino 3.4 index as one that has both 

360 maximum correlation with the core ENSO phenomenon and strongest influence on remote 

361 teleconnection events. It is the area-averaged SSTA over the region bounded by 5°N-5°S and 

362 120°W-170°W, straddling the Nino 3 and Nino 4 regions. Monthly data for the Nino 3.4 index was 

363 obtained from the International Research Institute (IRI) on Climate and Society Data Library 
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364 (http://iridl.ldeo.columbia.edU/SOURCES/.Indices/.nino/.EXTENDED/). Figure 3 gives a plot of the 

365 Nino 3.4 index. 

366 4.2.3 Multi-ENSO Index [MEI] 

367 MEI is not based solely on SST but is a multivariate index based on six variables recorded over the 

368 tropical Pacific and published in the Comprehensive Ocean-Atmosphere Data Set (COADS): sea- 

369 level pressure, zonal and meridional components of the surface wind, sea surface temperature, 

370 surface air temperature, and total cloudiness fraction of the sky. It is the first unrotated principal 

371 component of all the six observed fields combined [ Wolter and Timlin, 1993], and is analyzed 

372 separately for twelve sliding bi-monthly seasons, which removes most intra-seasonal noise. For 

373 correlation studies, it is advised that the MEI values for month (/ - 1) and month (/) be treated as for 

374 month (/). Monthly data for MEI for 1950 to 2008 was obtained from NOAA’s Earth System 

375 Research Laboratory (http://www.esrl.noaa.gOv/psd//people/klaus.wolter/MEEmei.html). Figure 3 

376 gives a plot of the MEI. 

377 4.2.4 El Nino Modoki [EMI] 

378 EMI is available from the Japan Agency for Marine-Earth Science and Technology 

379 (http://www.jamstec.go.jp/frcgc/research/dl/iod/modoki_home.html.en). It is defined as 

EMI = [SSTA\ c - O.S^STd],, - 0.5[AST4] w . (1) 

380 where [SSTA] is the area-averaged SST anomaly for the following regions: C (central, 165°E- 

381 140°W, 10°S-10°N), E (eastern, 110°W-70°W, 15°S-5°N), and W (western, 125°E-145°E, 10°S- 

382 20°N) [ Ashok et a!., 2007]. The time series for EMI is given in Figure 3. 
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383 Figure 3: Time series of PDO, Nino 3.4, MEI, and EMI indices. The PDO, MEI, and EMI 

384 series are overlain with a running, centered 13-month average to highlight yearly variations. 

385 The normalized monthly Nino 3.4 index series is overlain with a centered 3-month running 

386 mean and ± 0.5°C thresholds (see section 5.3.1). 

387 4.3 Forcing Data and Parameters for LSM 

388 The North American Land Data Assimilation System - Phase 2 (NLDAS-2) forcing data was used 

389 to run the Noah LSM. NLDAS-2 has a 1/8° latitude/longitude resolution over a domain covering 

390 CONUS, part of Canada and Mexico (125°W-67°W, 25°N-53°N), thus allowing the modeling of 

391 both the US and Mexican portion of RG. It incorporates both measured and modeled data from 

392 multiple sources: gauged precipitation measurements, satellite data, radar precipitation 

393 measurements, and output from numerical prediction models. NLDAS has been run retrospectively 

394 starting in January 1979, and provide hourly measurements in near real-time. The dataset include u 

395 and v wind components at 10 m above the surface, air temperature and specific humidity at 2 m 

396 above the surface, surface pressure, surface downward longwave and shortwave radiation, and total 

397 and convective fraction of precipitation, convective available potential energy, and potential 

398 evaporation. 

399 The precipitation field is from NCEP Climate Prediction Center (CPC) reprocessed daily gauged 

400 analyses that have been subjected to orographic adjustment based on the Parameter-elevation 

401 Regressions on Independent Slopes Model [PRISM; Daly et al., 1994] climatology and interpolated 

402 to the 1/8° NLDAS grid, and temporally disaggregated to hourly timescale using either the NWS 

403 Doppler radar-based (WSR-88D) precipitation, which has a 4 km spatial coverage, or the 8 km 

404 NOAA CPC Morphing Technique (CMORPH) hourly precipitation analyses [Joyce et al., 2004], 

405 This product uses hourly multiagency gauge data for bias correction and has been mosaicked over 

406 CONUS by NCEP/EMC [Baldwin and Mitchell, 1997]. The radar network is limited at the US 

407 borders with Canada and Mexico, and around 13% of CONUS are not covered. These gaps are first 
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408 filled with nearest neighbor mosaicked data from within a 2° radius, and if the latter is not available, 

409 CMORPH data are used instead. In Mexico, which is outside the radar covering range, CMORPH 

410 data is used. CMORPH is not available before 2002, and CCP hourly data is used and if it is not 

411 available, the North American Regional Reanalysis [NARR; Mesinger et al., 2006] data is used 

412 instead. 

413 The NLDAS dataset has been extensively compared, tested, and validated for snow cover and snow 

414 water equivalent [Pan et al., 2003; Sheffield et al., 2003], soil moisture [ Schaake et al., 2004], and 

415 streamflow and water balance [ Lohmann et al., 2004]. It has also been evaluated against the 

416 Atmospheric Measurement Program/cloud and radiation testbed, the Surface Radiation observation 

417 data, and the Oklahoma Mesonet, which is a quality controlled dataset from a dense network of over 

418 100 meteorological stations with meteorological measurements taken every 5 minutes [Luo et al., 

419 2003; Robock et al., 2003]. 

420 Noah also requires a set of parameters defining soil, vegetation, and topography for each grid. We 

421 use Zobler’s assessment of Food and Agriculture Organization (FAO) Soil Units [Zobler, 1986] 

422 which gives sand, silt, and clay fractions. The land cover is from the University of Maryland’s 

423 (UMD) 1 km spatial resolution global land cover product [Hansen et al, 2000]. It contains 14 

424 different classes (11 vegetation types, bare ground, urban/built up area, and water). The model also 

425 requires information on the quarterly and maximum albedo, monthly greenness fraction, and bottom 

426 temperature without elevation corrections. 
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427 5 RESULTS AND DISCUSSION 

428 5.1 Hydroclimatology of the Basin 

429 5.1.1 Spatial Variation 

430 RG trends across different climatic zones - alpine in southern Colorado and northern New Mexico, 

431 desert in southern New Mexico and west Texas, humid continental in east Texas, and humid sub- 

432 tropical in south Texas and north-eastern Mexico \Dahm et al., 2005] - making it an interesting 

433 study basin from both a hydrological and ecological perspective. Average annual precipitation varies 

434 from northwest to southeast across the basin, with a minimum of 187 mm at Manassa, in the San 

435 Juan Mountains, and a maximum of 698 mm at Port Isabel at the mouth of the basin. In the Rio 

436 Conchos sub-basin, precipitation varies from southwest to northeast with a maximum annual mean 

437 of 781 mm at El Vergel, in the Sierra Madre Occidental, and a minimum of around 290 mm around 

438 the mouth as it discharges into RG. Figure 4 shows the isohyet of the annual mean precipitation and 

439 the coefficient of variation (C v ) of monthly precipitation across the basin. C v is a statistical 

440 measure of variability, where a C v < 1 indicates less variation, while a C v > 1 indicates high 

441 variability. 

442 Figure 4: (a) Isohyet of annual mean precipitation and coefficient of variation of monthly 

443 precipitation across the Rio Grande basin, and (b) mean monthly precipitation for each region 

444 The basin exhibits wide disparity in precipitation regime, with C v ranging between 0.7 and 1.7. The 

445 Upper RG region has low variability. It receives around 20 mm of precipitation every month except 

446 for JAS with August being the wettest month (Figure 4(b)). The middle portion of the basin exhibits 

447 high variability. It has a unimodal precipitation regime, typical of the southwestern US where, for 

448 most of the year the average precipitation is below 20 mm except for May/June to 

449 September/October when the North American monsoon (NAM) brings most of the yearly rainfall. 
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450 Rio Conchos follows a similar pattern. NAM is a regional-scale circulation that develops over 

451 southwest North America, bringing substantial rainfall to this otherwise arid region. It is associated 

452 with a subtropical ridge shifting poleward during the summer months over the northwestern 

453 Mexican plateau and southwestern US. As evidenced in the mean monthly precipitation distribution 

454 (Figure 4(b)) of the Rio Conchos basin, the Upper-Middle and Middle-Middle RG, NAM starts to 

455 develop in late May to early June in southern Mexico quickly spreading along the western slopes of 

456 the Sierra Madre Occidental and into New Mexico and the western edge of Texas in early July and 

457 into southwestern US in the middle of July and decays in September/October [Adams and Comrie, 

458 1997; Higgins et a/., 1997]. The strength and path of the subtropical ridge is influenced by ENSO 

459 conditions. El Nino (La Nina) influences NAM by causing a weaker (stronger) southward 

460 (northward) displaced monsoon ridge [ Castro et al., 2001], The lower part of the basin has a slightly 

461 bimodal precipitation pattern with May and June and September and October as the two wettest 

462 periods. Atlantic hurricanes bring copious amounts of rainfall in a very short period of time, often 

463 resulting in major flooding in this region. A strong relationship between ENSO and the frequency of 

464 Atlantic hurricanes has been noted, with El Nino (La Nina) favoring a decrease (increase) in activity 

465 [Pielke and Landsea, 1999]. 

466 5.1.2 Temporal Variation 

467 Figure 5 shows the time series of the average monthly precipitation from all stations within each 

468 region, smoothed with a centered 13-month moving average window to accentuate yearly variations. 

469 Precipitation from Mexican stations in Lower RG is shown in a separate plot. The top panel of the 

470 figure shows the duration of PDO and ENSO in their respective phases. Yearly standardized 

471 anomaly, which is the anomaly for a particular year divided by the standard deviation of the series, 

472 was calculated for each station. A positive (negative) standardized anomaly indicates higher (lower) 

473 than average precipitation. The patterns give a visual picture of the time, duration, and spatial extent 
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474 of deficits and excesses in yearly precipitation across the basin. The standardized anomaly for the 

475 stations in the US, Mexico, and Rio Conchos are shown separately in Figure 6. 

476 Figure 5: Smoothed series of average monthly gauged precipitation (mm/month) for each 

477 region. The broken horizontal line is the long-term average. The scales of the ordinates are 

478 arbitrary to illustrate variability. The top panel shows the phases of PDO and ENSO. 


479 Figure 6: Time series of standardized anomaly for stations in the US, Mexico, and Rio 

480 Conchos (Stations are latitudinally arranged with numbers corresponding to Figure 2) 

481 Precipitation across the basin varies both spatially (Figure 4) and temporally. The standardized 

482 anomaly plot shows that precipitation is generally around the long term mean interspersed with 

483 drought spells, and a few extremely wet years, often spreading across the whole watershed. There 

484 have been more periods of rainfall deficits than excesses between 1935 and 2008. The most critical 

485 drought event started in 1951 and lasted up to 1956. In the semi-arid upper half portion of the basin, 

486 the drought started earlier, in 1942, right after an exceptionally wet year (1941). This event, 

487 commonly tenned as the drought of the 1950s, affected a large extent of the contenninous US and is 

488 the most severe drought on record for the watershed. It coincides with the cold phases of both ENSO 

489 and PDO. ENSO actually oscillated mostly between La Nina and neutral conditions between May 

490 1942 and February 1957 (Figure 3), the longest period in which the index remained in these phases. 

491 PDO shifted into its cold phase in November 1947 and remained mainly in this phase for up to 

492 February 1957. Another significant period of rainfall deficit started in 1962 and lasted up to 1965. It 

493 is associated with a drought that affected most of the northeast US [Namais, 1966; 1967] and 

494 coincides with cold PDO and mostly neutral ENSO with La Nina conditions between May 1964 and 

495 February 1965. 
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496 The 1976-1988 and 1991-1994 periods were that of nonnal and above normal precipitation across 

497 the watershed, with a few dry patches. Meteorological conditions prevailing were influenced by 

498 warm PDO (which started in July 1976) reinforcing multi-year El Ninos (see Table 1). Year 1989, 

499 however, was a dry year across the whole basin, sandwiched between two wet periods. This lack of 

500 precipitation was influenced by negative ENSO (1988-1989 La Nina) and a weak PDO oscillating 

501 mostly between ±0.5°C. Given that the 1988-1989 La Nina is one of the strongest (see Table 2), and 

502 is not reinforced (dampened) by a negative (positive) PDO, its influence on water availability is 

503 worth investigating further. 

504 The 2002-2004 drought is the second longest and most severe the region has experienced in recent 

505 record. It started earlier in the lower central arid portion and then propagated across the basin. 

506 Interestingly, this event did not coincide with La Nina conditions, but rather with neutral ENSO 

507 shifting into El Nino conditions, while PDO was in its cold phase until mid-2002 and shifted into its 

508 warm phase thereafter. Such dry conditions are not uncommon in the southwest as neutral ENSO 

509 with cold PDO can result in nearly as dry conditions as La Nina with cold PDO [ Goodrich , 2004; 

510 Quiring and Goodrich, 2008], The impact of such combination of the two indices on water 

511 availability will be further examined. 

512 5.2 Climate Teleconnection with Precipitation 

513 Precipitation is the main detenninant of water availability. To assess the relationship and spatial 

514 influence of climate variability on the precipitation in RG, climate indices were correlated with the 

515 average seasonal precipitation for each rainfall station. Kriging plots [ Delhomme , 1978] for 

516 correlations with PDO (Figure 7), and ENSO indices (Figure 8) were constructed for winter (DJF), 

517 spring (MAM), summer (JJA), and fall (SON). The Pearson correlation coefficient, p , forms the 

518 basis of the statistical test of independence. The null hypothesis is that the seasonal average 
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519 


precipitation is iid nonnal random variable not dependent on the indices. The magnitude and sign of 

520 the correlation coefficient are thus indicative of the existence, strength, and nature of any association 

521 [Redmond and Koch, 1991]. The correlation significance (p-value) for each station was also 

522 computed and is given (inset maps) along with the kriging plots. The /;- value is the probability, if 

523 the test statistics is distributed as assumed in the null hypothesis, of observing a test statistic as 

524 extreme, or more, than the one actually observed, i.e., how unlikely it may be due to chance. The p- 

525 values are stratified as follows: less than or equal to 0.01%, 0.1%, 1%, 1-5%, and greater than 5%. 

526 5.2.1 PDO 

521 The PDO series between 1935 and 2008 has almost equal number of cold and wann years, thus 

528 prevents any bias that may result from the dominance of either phase over the other [ Kurtzman and 

529 Scanlon, 2007]. The correlation structure is different for each season and varies considerably across 

530 the basin (Figure 7). The highest correlations are in winter and spring. In winter, a statistically 

531 significant positive correlation, ranging between 0.3 and 0.5, is observed in the lower RG region, 

532 especially in Mexico, at the mouth of the basin. This correlation structure shifts diametrically in the 

533 spring, with the northern semi-arid regions exhibiting higher correlation, while those in the lower 

534 half, from the Rio Conchos downwards, are statistically insignificant. The stations in the Rio 

535 Conchos do not show any statistically significant correlation with PDO, except for one station in 

536 winter, at the mouth of this sub-basin. Summer precipitation, across the whole basin, exhibits no 

537 significant correlation with PDO, while negligible correlation is observed in some sections of the 

538 basin in the fall. Hence it is noted that despite relatively low correlation, PDO does have an 

539 influence on the winter and spring conditions in the southern and northern parts of the basin 

540 respectively. Therefore knowledge of the state of PDO is essential for water resources planning. 
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541 5.2.2 ENSO - Nino 3.4, MEI, and EMI 

542 The seasonal kriging plots for Nino 3.4 is given in Figure 8 and those for MEI and EMI are 

543 available as auxiliary material to this paper. All three indices exhibit the same general seasonal 

544 correlation structure, but overall the correlation of EMI with precipitation is the weakest. 

545 Just like for PDO the correlation structure is different for each season and varies across the basin. A 

546 positive, statistically significant correlation, ranging between 0.3 and 0.7, was observed for winter 

547 and spring across the whole basin, except for the Upper RG region which exhibits a negative 

548 correlation with ENSO. This negative correlation at the head of basin, in the Rockies, is consistent 

549 with the findings of Smith and O'Brien [2001] and Patten et al. [2003] on snowfall pattern in the 

550 US. In spring the highest correlation region is the Upper-Middle RG. From a water management 

551 perspective this finding is important as higher snowfall at the headwaters, during La Nina 

552 conditions, may offset reduced precipitation in New Mexico. In summer and fall a negative 

553 correlation, even though not statistically significant, is observed, especially in the lower half of the 

554 basin. 

555 The correlation between EMI and precipitation is not statistically representative in the basin and 

556 may therefore not be useful for water management purposes. The correlation structures of Nino 3.4 

557 and MEI are similar; therefore despite MEI’s appeal as a more inclusive ENSO index, we shall 

558 adhere to NOAA’s operational definition of El Nino and La Nina, as discussed in section 5.3.1 

559 below, for the remainder of the paper. 

560 Figure 7: Plots of seasonal correlation coefficients between PDO and precipitation anomaly, 

561 Inset gives the /7-values for the regression coefficients. 
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Figure 8: Plots of seasonal correlation coefficients between Nino 3.4 and precipitation anomaly 
(Similar plots for MEI and EMI are available as auxiliary material). Inset gives the /7-values 

for the regression coefficients. 


5.3 ENSO Events 

5.3.1 Definition of El Nino and La Nina 

The definition for El Nino has evolved over time; different investigations use different indices and 
criteria, thus producing dissimilar lists of events. For long, there was no specific definition for La 
Nina despite an ongoing debate in the scientific community; they were defined in the context of the 
El Nino phenomenon [ O'Brien , 2002]. Trenberth [1997] analyzed El Nino conditions between 1950 
and 1997 using both Nino 3 and Nino 3.4 indices, relative to a base period climatology of 1950- 
1979, and suggested that an ENSO event is deemed to occur when the Nino 3.4 index is above (or 
below) ±0.4°C for at least six months. In 2003 the National Oceanic and Atmospheric 
Administration (NOAA) issued an official operational definition for El Nino and La Nina. The 
definition was endorsed by both Canada and Mexico in 2005 (North American countries reach 
consensus on El Nino definition available at http://www.noaanews.noaa.gov/stories2005/s2394.htm) 
and has been adopted by the World Meteorological Organization Region IV [Larkin and Harrison, 
2005]. Accordingly, an El Nino (La Nina) is defined as: 

A phenomenon in the equatorial Pacific Ocean characterized by a positive (negative) sea 
surface temperature departure from nonnal (for the 1971-2000 base period) in the Nino 
3.4 region greater than or equal in magnitude to 0.5 degrees C, averaged over three 
consecutive months. 

Neutral condition is when the index is between ±0.5°C. Figure 3 shows a plot of the Nino 3.4 index 
for 1935 to 2008, overlain with a centered 3-month moving average, which smoothens out variations 
in SST not associated with ENSO, and the ±0.5°C thresholds. The base period climatology for the 
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586 IRI dataset is 1951-1980 [Kaplan et al., 1998] and was adjusted to 1971-2000 to satisfy the above 

587 definition. During that period, El Nino occurred 21% of the time and La Nina 27% of the time, i.e. 

588 the Pacific was either active in one of the two conditions or neutral about half the time. 

589 5.3.2 Comparison ofENSO Events (1 979 - 2008) 

590 Table 1 gives the start, end, and duration of El Nino and La Nina events occurring since 1979. The 

591 braces indicate coupled events where the average SSTA remained in one phase. NOAA’s definition, 

592 unlike Trenberth ’s [1997], does not specify a minimum period over which the average index has to 

593 be above (or below) the threshold to be deemed a significant event, thus smaller periods of one to 

594 three months are also included in the list. We note that between 1979 and 2008, El Ninos generally 

595 lasted longer than La Ninas, and there were more months above the +0.5°C threshold, resulting in El 

596 Ninos 26.4% of the time, as opposed to 21.1% for La Ninas. A higher frequency of El Ninos since 

597 1976 has been associated with decadal changes in the climate in the Pacific Ocean accentuated by 

598 recent climatic changes [Kumar et al., 1994; Trenberth, 1990]. The two often called major El Ninos 

599 (1991-1995 and 2002-2005) and major La Ninas (1983-1986 and 1998-2001) occurred between 

600 1979 and 2008. 

601 Table 1: El Nino and La Nina events between 1979-2008 following NOAA’s definition 

602 The onset of ENSO event is not consistent; El Nino events typically begin between May and 

603 September and La Nina events begin between July and October, except for the 1988 event which 

604 started earlier, in May, and the 2005 event evolved later, in December. The time interval between El 

605 Nino events ranged between 15 and 47 months with an average of 33 months when coupled events 

606 are counted as one. In the case of La Nina events, the interval was normally around 26 months, 

607 except when the two major El Nino events occurred, leading to a gap of 75 and 57 months 
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608 respectively. The time interval between phase shifts varies from as short as one month (1997-1998 

609 El Nino to 1998-2001 La Nina) to neutral conditions lingering for up to 23 months (1988-1989 La 

610 Nina to 1991-1995 El Nino). 

611 Not all El Ninos (La Ninas) are same; some events last for only a few months while others may 

612 persist for two or more calendar years. Figure 9 compares the recent strongest El Nino and La Nina 

613 events. El Ninos (La Ninas) tend to peak (trough) between November and February, except for the 

614 1986-1988 event which peaked earlier in September. The strongest El Nino has a higher Nino 3.4 

615 index value compared to the lowest index value for the strongest La Nina. Two El Ninos with 

616 maximum SSTA above +2°C were noted while La Ninas’ maxima were within -1 and -2°C. The 

617 reason for this lies mainly in the depth of the thermocline, which is shallower on the easternmost 

618 part of the Pacific basin, and it is therefore harder for the cold tongue to get colder, thus explaining 

619 why El Ninos tend to be stronger than La Ninas [ Neelin , 2010]. 

620 All these differences make it difficult to rank El Ninos (La Ninas). It is hard to define the criteria to 

621 be used for such ranking; should it be (i) the duration above (or below) the defined threshold, (ii) 

622 the maximum (or minimum) SST recorded, or (iii) the total duration of the event in the phase of 

623 interest [ Wolter and Timlin, 1998]? Glantz [1998] argues that a number of socio economic criteria, 

624 such as global spread of impacts, costs of devastation, or even public perception and media 

625 coverage, should also be considered. These criteria, however, are too complex and not totally 

626 objective and will thus not be considered in this study. Table 2 ranks the three strongest ENSO 

627 events, using the Nino 3.4 index, following Wolter and Timlin [1998], and intensity - a new metric 

628 which incorporates both SSTA and duration. ENSO intensity is the sum of monthly indices above 

629 (or below) the ±0.5°C threshold divided by the corresponding number of months. For combined 

630 events, the months where the index is within the neutral range are excluded. 
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631 Figure 9: Comparison of major El Nino and La Nina events between 1979 and 2008 

632 (transparent gray band shows ±0.5°C neutral range) 

633 Table 2: The three strongest ENSO events since 1979 based on different ranking criteria 

634 Duration and Maximum/Minimum SSTA 

635 The ranking reveals some interesting results regarding the individuality of ENSO events. In both El 

636 Ninos and La Ninas, the list of the three strongest events using duration above the threshold is very 

637 different from that using the maximum (or minimum), implying that events persisting for long 

638 periods do not necessarily exhibit large deviations from the mean. The 1991-1995 event, for 

639 example, started in January 1990 and persisted until July 1995 [ Trenberth and Hoar, 1996], and is 

640 longest El Nino on record, with a combined total of 28 months above the +0.5°C threshold, but does 

641 not show up within the list of events with maximum temperature or intensity. Allan and D'Arrigo 

642 [1999], examining palaeoclimatic records of ENSO infer that such events are not unusual. El Ninos 

643 persisting for three years or longer may occur around four to five times per century. On the other 

644 hand, instrumental records analyzed by Trenberth and Hoar [1996] indicate that this event is 

645 unlikely the result of natural decadal-timescale variation but is rather influenced by global wanning. 

646 The 1997-1998 El Nino, which developed earlier than the scientific community expected, is also 

647 noteworthy and has been dubbed “the climate event of the century” [ Changnon and Bell, 2000], 

648 Unlike the 1991-1995 event, it lasted only 13 months but was still one of the strongest on record 

649 with a peak SSTA of 2.56°C, closely rivaling the 1982-1983 event [Kiladis and Diaz, 1986]. By 

650 contrast the 1990-1995 event had a peak temperature of only 1.71°C occurring within the first third 

651 of the event. This super El Nino, just like the one in 1982-1983, had a dramatic impact on global 

652 weather variability causing the second worst drought in Australia and devastating floods across the 

653 western United States, with record precipitation in California. Western Pacific Islands, Central 
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654 America and Mexico experienced severe droughts. It ended abruptly in the middle of 1998 when a 

655 mass of cold water under-rode the wann surface waters causing one of the longest and strongest La 

656 Nina ever recorded. 

657 The 1998-2001 La Nina started in June 1998 and lingered till May 2001, with a combined total of 27 

658 months below the -0.5°C threshold. The SSTA reached a minimum of -1.50°C in December 1998 

659 and appeared to decay thereafter but resurged with a lower SSTA of -1.76°C in January 2000. The 

660 1988-1989 La Nina, the second longest since 1979, reached a lower SSTA of -2.03°C in November 

661 1988, and is one of the lowest on record. 

662 When the total duration the index remain in the positive or negative phase is considered, the same El 

663 Ninos are ranked, but a slight difference is noted in the order of La Ninas. The 1983-1986 event, 

664 even though has fewer months below the threshold, takes longer to decay, and hence supersedes the 

665 1988-1989 event. 

666 ENSO Intensity 

667 ENSO intensity shuffles the ranking of El Ninos based on maximum SSTA recorded. The 1997- 

668 1998 event, with an intensity of 1.325, surpasses the 1982-1983 event, even though the latter has the 

669 highest temperature among recent episodes. In the case of La Nina, the list based on intensity agrees 

670 with that based on temperature, but not with the one based on duration. The 2007-2008 La Nina, 

671 being a short event, does not feature in the list based on duration but has the second lowest 

672 temperature and intensity. It is also important to note that La Ninas have much higher intensities 

673 than El Ninos as they generally have shorter durations. 
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The effect of ENSO on water availability, based on the criteria discussed above, will be examined to 
ascertain if SSTA, duration, or intensity have different impacts. Note that in subsequent sections, 
ENSO events will be referred to as labeled in the legend of Figure 9. 

5 .4 Runoff and Water Availability 

Noah or LIS does not have a routing scheme; runoff is generated at 1/8° pixel level. Further, since 
the basin is divided into sub-regions rather than sub-basins, and because RG has a number of 
endorheic sub-basins, area-averaged runoff (AAR) is a better representation of regional water 
availability [ Shukla and Wood, 2008], 


AAR = 



r : 


( 2 ) 


where r is the pixel runoff at the temporal scale of interest, a is the area of the corresponding pixel, 
and n is the number of pixel in the sub-region considered. 

Before assessing changes in water availability, we evaluated NLDAS-2 and runoff from Noah LSM 
against observations to assess the validity of our approach. 

5.4.1 Model Validation 
Precipitation 

The precipitation field in NLDAS-2 is derived from a number of sources whose coverage near the 
US borders and beyond is often limited. The product has been extensively validated over CONUS 
but has not received the same treatment beyond the US borders. Given that RG is a transboundary 
basin, extending into Mexico, and that the density and length of precipitation record in Mexico is 
low, it is important to verify the adequacy of NLDAS-2 over the Mexican portion of the basin. We 
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selected three of the six sub-regions within the basin, one in the US and two in Mexico, for 
validation. The Upper-Middle RG has a dense rain gauge network, with 31 stations; it is found 
almost entirely within the state of New Mexico and therefore benefits from the high quality products 
used in generating NLDAS-2. The two regions in Mexico are the Rio Conchos and the Lower RG. 
Rio Conchos, located entirely within Mexico, is an important sub-basin, but has only four viable 
stations aligned along the major axis of the catchment. The major portion of Lower RG is located 
within Mexico, and has 5 stations along the US border and 11 stations in Mexico. We therefore 
compared the NLDAS-2 precipitation field with the observations from the Mexican stations. 

The precipitation field in NLDAS does not agree very well with station observations at the hourly 
timescale, but as the data is aggregated over longer timescales the correlation increases [Luo et ah, 
2003]. This is because NLDAS-2 precipitation is generated from multiple sources and averaged over 
the domain of interest. We compared precipitation at the monthly scale. Figure 10 shows the scatter 
plots of the area-averaged monthly precipitation (similar to AAR) derived from the NLDAS-2 
precipitation field (m) against average monthly observed precipitation from gauges (O) within 
each sub-region. The data length differs, as shown by the total number of data points (N) in each 
plot. The relative bias, root mean square deviation (RMSD), and correlation coefficient ( p ) are 
given for each region, along with the identity (1:1) line. 

Relative Bias = ^ ® ( 3 ) 


RMSD = 


j Z><-A 2 

V N 
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710 Figure 10: Scatter plot of area-averaged monthly NLDAS-2 precipitation and observed 

711 precipitation for the (a) Upper-Middle RG, (b) Rio Conchos, and (c) Lower RG region 

712 (Mexican stations only) 

713 All three regions have very high p values implying that the precipitation field in NLDAS-2, at the 

714 monthly timescale, over the entire RG basin is in good agreement with observations. A small 

715 consistent bias is noted towards observations, similar to the findings reported by Luo et al. [2003]. 

716 The RMSD values for the Rio Conchos and Lower RG are greater than that for the Upper-Middle 

717 RG by a factor of three. This can be attributed to the fact that the Upper-Middle RG has a dense 

718 network unifonnly spread over the whole region, while in Mexico the density and spread is limited. 

719 Also, most stations in Mexico did not have a continuous series but have been reconstructed from 

720 observations available from adjoining sites. 

721 Runoff 

722 In order to assess the representativeness of the modeled runoff we compared AAR of the Rio 

723 Conchos, which is a closed sub-basin, unlike the other sub-regions, to measured values at Ojinaga, 

724 which is located at the confluence of the Rio Conchos with RG. Historical mean daily discharge data 

725 for RG at Ojinaga is available from the International Boundary and Water Commission (IBWC) 

726 (http://www.ibwc.state.gov/). This dataset, as can be expected, incorporates land-use-land-cover 

727 changes, the effect of dams, diversions, and other infrastructural changes. A dataset of naturalized 

728 monthly flow, extending up to 2000, for the Rio Conchos at Ojinaga is also available \Sandoval- 

729 Solis et al, 2010]. The naturalization process utilized streamflow recorded at gaging stations and 

730 adjusted to remove the effect of reservoir storage and evaporation, water supply diversions, and 

731 return flows from surface and groundwater, so that the resulting series is as close to flow unimpaired 

732 by engineering infrastructures. However the influence of land-use-land-cover changes, infiltration, 

733 surface storage-flow, subsurface storage-flow, and evapotranspiration cannot be adequately 
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734 accounted for \Wurbs, 2006]. Noah accounts for evaporation, infiltration, and other hydrological 

735 parameters; but land-use-land-cover was kept constant to minimize extraneous noise that may 

736 interfere with the climate teleconnection signals. Figure 11 (a) gives the plot of AAR, naturalized, 

737 and observed time series of monthly streamflow in the Rio Conchos. Prior to any inference, it is 

738 important to point out that none of these three series are accurate representation of flow in the basin; 

739 each have different intrinsic limitations. Figure 11(b) gives a 12 -month- sliding- window correlation 

740 plot for the period January 1998 to December 2000, which includes the period for which Noah was 

741 validated by Lohmann et al. [2004], thereby giving a benchmark for comparison. Each value in the 

742 graph is the Pearson correlation coefficient for 12 consecutive months starting from the month at 

743 which it is plotted. This process allows us to compare the consistency in the correlation across 

744 different seasons. 

745 Figure 11: (a) Modeled, naturalized, and measured streamflow in the Rio Conchos sub-basin, 

746 (b) sliding-window correlation 

747 It can be noted that Noah faithfully captures the monthly variations in runoff in the basin. 

748 Surprisingly, however, the correlation between AAR and the measured flow is consistently higher 

749 than with naturalized flow. This may be a function of the naturalization procedure, where 

750 evapotranspiration and other hydrologic process cannot be amply determined given their 

751 complexity. Another notable feature of the running correlation is the periodic variation in the 

752 correlation values. The correlation is generally between 0.85 and 0.95 except for certain short 

753 intervals, which may be attributed to reservoir operation and diversions affecting the recorded flow. 

754 Lohmann et al. [2004] reported a correlation of 0.954 for the Nehalem River in Oregon which is a 

755 much smaller watershed (1,727 km") compared to the Rio Conchos (64,000 km"). In the larger snow 

756 dominated Wind River in Wyoming (4,897 km"), a very low correlation of 0.117 was obtained 
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757 which was due to a difference in the timing of snowmelt in Noah. The highest correlation across the 

758 series is generally noted in the fall and winter and this can be explained by multiple factors: the 

759 timing of snowmelt [Pan et al., 2003; Sheffield et al., 2003], the influence of the North American 

760 Monsoon [Hogue et al., 2005], and the influence of engineering infrastructures, such as dams, return 

761 flows, and unmetered diversions for agricultural purposes within the basin. Given that summer and 

762 early fall are the wettest months for the region (Figure 4), streamflow may initially stored in the 

763 reservoirs and any excess thereafter is released once the reservoirs have reached capacity and hence 

764 recorded at the downstream stream gauge as values closest to the natural variations. 

765 5.4.2 Temporal Pattern - Climate Indices, Precipitation, and Water Availability 

766 Continuous wavelet transform allows the study of the temporal structure of precipitation and runoff 

767 across the basin and make inferences on the influence of climate variability patterns. The continuous 

768 wavelet power spectrum of AAR from each section of the basin is given in Figure 12(b) and can be 

769 compared against those of the four climate indices considered in this study (Figure 12(a)). The time 

770 spans for the plots of climate indices are the same as for runoff, i.e., 1979 to 2008. We used the 

771 MATLAB" based software package developed by Grinsted et al. [2004] (available from 

772 http://www.pol.ac.uk/home/research/waveletcoherence/) for generating the wavelet plots. The 

773 statistical significance of the peaks in the wavelet spectrum was tested using Monte Carlo methods 

774 against a lag-1 autoregressive red noise background; peaks with greater than 95% confidence are 

775 designated by thick black contours in the figures. However, regions out of these 95% confidence 

776 level areas should not be construed as the product of noise only. Natural processes are also present 

777 in these regions, albeit having a lower bearing on the power spectrum, and infonnation on the 

778 influence of climate teleconnections can still be garnered from them [Anctil and Coulibaly, 2004]. 
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779 Figure 12: Continuous wavelet power spectrum of (a) climate indices and (b) AAR for each 

780 region. Period is in days. The thick black contours designate the 5% significance level against 

781 red noise and the thin black line demarcate the cone of influence beyond which, shown in a 

782 lighter shade, the image may be distorted since the data is not padded at the edges. 

783 PDO has a cycle of about 20 to 30 years, which is not apparent in this wavelet power spectrum, as it 

784 is limited to 1979 to 2008 only. During this time span, the index was mostly positive until 1995 after 

785 which it oscillated with a period of three and a half years in each phase (Figure 3). This cycle is 

786 visible in the continuous wavelet power spectrum. ENSO has a much shorter wavelength, with a 

787 recurrence pattern of 3 to 6 years and every event normally lasts for around a year. Significant 

788 power in this band is observed throughout the entire record in the wavelet power spectrum of both 

789 Nino 3.4 and MEI, along with weaker significant power in the 1.5 to 2 year band associated with 

790 secondary variations in the indices while in a particular phase. Nino 3.4 and MEI have very similar 

791 power spectrum as they are highly correlated (p = 0.92). The plot for EMI is different from the latter 

792 indices as it has a different mode and evolves with a different frequency [Ashok et al., 2007]. 

793 Variations at both ENSO and PDO frequencies are apparent in the wavelet power spectrum of AAR. 

794 Significant powers in the smaller period (higher frequency) band are the result of seasonal 

795 variations. It is worth pointing that the continuous wavelet power spectra of both the averaged 

796 gauged precipitation and the first principal component of gauged precipitation for each sub-region 

797 (not shown here) exhibited variations similar to ENSO and PDO frequencies, further supporting the 

798 correlations established in section 5.2. The fact that the influence of climate variability was more 

799 apparent in precipitation than in streamflow is because runoff is not a first order response of 

800 precipitation but is filtered by the watershed characteristics [ Legates et al., 2005], Looking closer at 

801 each spectrum we note that in the Upper RG region significant power at the 5% significance level is 

802 exhibited within the high frequency band, coincident with the pattern exhibited by the PDO 

803 spectrum but ENSO related patterns were not significant at the 5% level. In the Upper-Middle RG 
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804 region, however, variations at the ENSO frequencies are clearly evident in most part of the spectrum 

805 except for the period 1995 to 2003. In the Middle-Middle RG region, the spectrum exhibited 

806 variations in the ENSO frequencies even though not statistically significant throughout. The 

807 spectrum for the Lower-Middle RG region was devoid of statistically significant powers as can be 

808 expected, because the correlation of ENSO indices with precipitation is consistently lower than in its 

809 two adjoining regions (Figure 8). In the Lower RG region variations at ENSO frequencies is visible 

810 in the second half of the series. It should be noted that this region is close to the Gulf of Mexico and 

811 subjected to oceanic influences and hurricanes whose effects are embedded in the model outputs. In 

812 the Rio Conchos, variations at ENSO frequencies are evident in the first half of the series but 

813 completely absent in the second part when the basin was subjected to an exceptional drought. 

814 5.5 Effects of Large-Scale Climate Indices on Water Availability 

815 After having established the spatial correlation of precipitation with large scale climate indices and 

816 investigated the sensitivity of runoff to PDO and ENSO using continuous wavelet transfonn, we 

817 now discuss the changes in water availability with respect to major El Nino and La Nina events. We 

818 also determine if there is any lag in runoff relative to ENSO. 

819 5.5.1 Lag Correlation 

820 To investigate the lag between individual ENSO events and runoff in RG we considered an 18 

821 months long series encompassing each event. Most events peaked in early winter; therefore the 

822 series considered extended from March before to August after the peak. The 1986-1988 El Nino 

823 peaked earlier in September and thus the range considered was from December before to May after 

824 the peak. Table 3 gives monthly lag correlation coefficients of AAR of the whole RG basin relative 

825 to the Nino 3.4 index. The lag is defined by the month having maximum statistically significant 

826 correlation. A positive lag indicates that the index led runoff. Lags ranging between 0 to 3 months 
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827 were observed for AAR of the whole basin with respect to most El Ninos and La Ninas. It is 

828 interesting to note that for some events a negative correlation was observed despite the fact that 

829 there is an overall dominant positive correlation between ENSO and runoff in the basin. 

830 Lag correlation was also computed for each sub-region and the consolidated result of the statistically 

831 significant correlations is shown in Table 4. For the 1998-2000 La Nina, which is bimodal, the lags 

832 shown are with respect to the second dip as it has a lower SSTA. In the Upper RG the lag is 3 to 4 

833 months. In the Middle-Middle RG a consistent lag of 1 to 4 months is observed relative to most El 

834 Nino events and in the Rio Conchos a lag of 0 to 5 months is observed relative to three of the four 

835 La Ninas considered. This result is consistent with the findings of Chen and Kumar [2002] who used 

836 a large-area basin-scale LSM to investigate the relationship among terrestrial hydrologic processes 

837 with ENSO over North America and Kumar and Hoerling [2003] who compared and confirmed the 

838 observed lag in zonal mean tropical thermal anomalies with respect to east Pacific SST using an 

839 atmospheric GCM. 

840 Table 3: Lag correlation of Nino 3.4 index with Noah runoff for the whole of the basin 

841 Table 4: Lag correlation for each section of the basin 

842 5.5.2 Changes in Water Availability 

843 In order to assess the effects of ENSO events on water availability, the seasonal percentage change 

844 with respect to the long tenn average in AAR, for each sub-region and for the whole basin was 

845 computed (Figure 13). ENSO events typically peak around November and since the lag between 

846 ENSO and runoff was found to be generally between 0 month to a season, water availability in 

847 winter and spring are expected to be most influenced. Note, however, that the 1986-1988 El Nino 

848 and 1998-2000 La Nina do not fit the general trend. The 1986-1988 El Nino lingered for a year 
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849 before reaching its peak temperature and the 1998-2000 La Nina was bi-modal, with two distinct 

850 troughs, over a 27-month period. 

851 The first thing we note is that even though there is a general tendency for an increase (decrease) in 

852 runoff during El Ninos (La Ninas), some events actually caused a decrease (increase) in water 

853 availability. 

854 Figure 13: Percentage change in water availability during (a) El Nino and (b) La Nina events 

855 PDO generally enhances ENSO events [Kurtzman and Scanlon, 2007], therefore we also consider 

856 the phases of PDO for corresponding major El Nino and La Nina events in our discussion. Table 5 

857 gives the coincident phases of PDO with respect to the ENSO events considered. We note that some 

858 El Ninos (La Ninas) were strengthened by positive (negative) PDO, while others coincided with 

859 weak or transitioning PDO. 

860 Table 5: PDO phase for major ENSO events 

861 El Nino and PDO 

862 The 1986-1988 El Nino, even though was the third strongest event since 1979 (Table 2), brought the 

863 highest percentage increase in runoff for the whole basin (196%), with the upper half of the basin 

864 gaining between 190 to 350% more runoff in winter. The lower half, including Rio Conchos 

865 experienced a more modest increase of around 60%. The same pattern, but with lower percentages, 

866 persisted in spring and summer. The event was enhanced by a strong positive PDO - both events 

867 evolved synchronously and peaked in August 1987 (Figure 3). 
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868 The 1991-1993 El Nino coincided with a PDO transitioning into its positive phase. It triggered the 

869 same pattern in runoff but with slightly lower percentages (average of 147% for the whole basin). 

870 The Rio Conchos bcnefitted from the highest increase in runoff in spring during that event, 

871 compared to all other El Ninos. This event had the longest duration above the ±0.5°C threshold 

872 among the El Ninos considered, but was not within the top three in terms of maximum temperature 

873 recorded. It had three peaks; the two following the first one, however, have lower SSTAs. The 

874 second peak coincided with positive PDO and the third one coincided with a temporary shift in the 

875 phase of PDO to negative. Interestingly, the percentage change in AAR was negative across the 

876 basin during both the second and third peaks. 

877 The 1982-1984, 1997-1999, and 2002-2004 El Ninos did not generate higher runoff in the basin, but 

878 rather a decrease in water availability. AAR was lower by 50% over the whole basin and over 90% 

879 in the Rio Conchos sub-basin. The 1982-1984 and 1997-1999 events were two of the strongest El 

880 Ninos on record, based on ENSO intensity and maximum temperature recorded. They both 

881 coincided with weak PDO; the 1982-1984 event starts when PDO was in its positive phase but 

882 peaked when the latter was almost neutral, oscillating between phases. The 1997-1999 event 

883 coincided with a PDO decaying from its positive to negative phase. 

884 La Nina and PDO 

885 AAR in the basin was normally lower than the long tenn average during La Nina winters and 

886 springs, except for the 1984-1986 event. The three events (1988-1990, 1998-2000, and 2007-2008) 

887 that caused a decrease in water availability occurred when PDO was either in its negative phase or 

888 transitioning from positive to negative. Given that the 1998-2000 event is the longest La Nina since 

889 1979 (27 months) and is bimodal, with the second dip, occurring in January 2000, having a lower 

890 SSTA and coinciding with negative PDO, we also computed the change in water availability 
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891 


following the second dip (not shown in Figure 13). A decrease in AAR was noted across all sub- 

892 regions, except for winter in Upper RG (increase of 6%) which exhibits negative correlation with 

893 ENSO. The total decrease in AAR for the whole basin was nearly 60% for year 2000. 

894 The 1984-1986 La Nina coincided with a positive PDO and was the only event that caused an 

895 increase in water availability. An exceptionally large increase is noted in the Rio Conchos in both 

896 winter and spring, which translated in an increase for the whole basin. These results are consistent 

897 with the analysis of multiyear droughts in the past three centuries by Cole et al. [2002] who showed 

898 that persistent negative PDO enhances the impact of La Nina related droughts, while oscillating 

899 PDO produced moderate and/or localized droughts, and a positive PDO would suppress drought 

900 despite persistent La Nina conditions. 

901 Based on the above observations, we note that PDO has an important influence on water availability 

902 in the basin. A positive PDO enhances the effect of El Nino and dampens the negative effect of La 

903 Nina. When PDO is in a neutral/transition phase La Nina dominates climatic conditions and reduces 

904 water availability. El Ninos lingering for long periods have more influence on water availability than 

905 short duration high intensity events. Finally it is interesting to note that the percentage increase 

906 during El Ninos significantly offsets the drought-causing effect of La Ninas. This finding should not 

907 be discounted in long-term water resources planning. 

908 6 CONCLUSIONS 

909 Local meteorological and hydrological variables, and hence water availability, are influenced by 

910 large-scale climate indices. In this study we investigated the influence of ENSO and PDO on the 

911 water availability in RG by first establishing the spatial and temporal variation of the correlation 

912 between climate indices and gauged precipitation across the basin and then detennining percentage 

913 changes in water availability as derived from an LSM instead of using streamflow which is 
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constantly impacted by activities in the basin masking climate influences. The following 

conclusions are drawn from this study: 

1. The correlation between PDO and three ENSO indices, namely Nino 3.4, MEI, and EMI, with 
gauged precipitation respectively shows that both ENSO and PDO have a strong influence on 
the winter and spring precipitation in the basin. The overall correlation is positive, except for 
the Upper RG section which includes the headwaters in the San Juan range in the Rocky 
Mountains in southern Colorado. Therefore, higher snowfall during La Nina conditions may 
help in maintaining flow in the river and offset precipitation reduction in arid/semi-arid New 
Mexico. 

2. The correlation between the Nino 3.4 and MEI indices with precipitation are similar since they 
are closely related. The temporal structure and influence of EMI is different and is not strongly 
correlated with precipitation in the basin. 

3. Additional infonnation can be garnered by examining the major El Nino and La Nina events 
by classifying them using four criteria (duration above defined threshold, maximum or 
minimum SSTA, duration in phase of interest, and intensity). ENSO events are not equivalent, 
some events have short duration but high intensity, while others lingers for several years with 
lower SSTA. 

4. Runoff across the basin was generated using the Noah LSM and AAR was used as a proxy for 
water availability. The basin was divided into six sub-regions for analysis purposes. 
Continuous wavelet power spectrum shows the extent of influence of ENSO and PDO on 
runoff. Variations at both ENSO and PDO frequencies are apparent in the wavelet power 
spectrum of AAR for each region. 
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5. The influence of individual ENSO events, five El Ninos and four La Ninas between 1979 and 
2008 and corresponding phases of PDO, on water availability in the basin was investigated. 
Lags ranging between 0 to 3 months were observed between runoff and ENSO events. A 
general increase (decrease) in runoff during El Ninos (La Ninas) was noted but some 
individual events actually caused a decrease (increase) in water availability. El Ninos lingering 
for long periods have more influence on water availability than short duration high intensity 
events. The upper-middle section of the basin records a higher increase in winter water 
availability during El Nino events (200-300%) while the lower half, including the Rio 
Conchos, experiences a more modest change. 

6. PDO has an important influence on water availability. A positive PDO enhances the effect of 
El Nino and dampens the negative effect of La Nina. When it is in its neutral/transition phase 
La Nina dominates climatic conditions and reduces water availability. 

7. The percentage increase during El Ninos significantly offset the decrease registered during La 
Ninas. This finding is important for water resources planning. 

The study extends the discussion between the influence of large-scale circulation indices and local 
meteorological and hydrological conditions by quantifying the seasonal percentage changes in water 
availability, which is more tangible information for water planning. Climate change may alter the 
frequency and intensity of ENSO events and may cause droughts that are more extreme and/or of 
longer duration than on record. The current results, while are not intended for prediction purposes, 
may help in the long-term sustainable water planning and management within the basin for both the 
United States and Mexico. Finally, the methodology adopted in this paper is not limited to the 
watershed scale but can be applied to larger continental scale to assess the need and effectiveness of 
interstate water transfers. 
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